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Mammalian genomes produce huge numbers of noncoding RNAs (ncRNAs). However, the functions of most ncRNAs are 
unclear, and novel techniques that can distinguish functional ncRNAs are needed. Studies of mRNAs have revealed that 
the half-life of each mRNA is closely related to its physiological function, raising the possibility that the RNA stability of 
an ncRNA reflects its function. In this study, we first determined the half-lives of 11,052 mRNAs and 1418 ncRNAs in HeLa 
Tet-off [TO) cells by developing a novel genome-wide method, which we named 5'-bromo-uridine immunoprecipitation 
chase-deep sequencing analysis [BRIC-seq). This method involved pulse-labeling endogenous RNAs with 5'-bromo-uridine 
and measuring the ongoing decrease in RNA levels over time using multifaceted deep sequencing. By analyzing the re- 
lationship between RNA half-lives and functional categories, we found that RNAs with a long half-life [f 1/2 > 4 h) contained 
a significant proportion of ncRNAs, as well as mRNAs involved in housekeeping functions, whereas RNAs with a short half- 
life [t\ /2 < 4 h) included known regulatory ncRNAs and regulatory mRNAs. The stabilities of a significant set of short-lived 
ncRNAs are regulated by external stimuli, such as retinoic acid treatment. In particular, we identified and characterized 
several novel long ncRNAs involved in cell proliferation from the group of short-lived ncRNAs. We designated this novel 
class of ncRNAs with a short half-life as Short-Lived noncoding Transcripts (SLiTs). We propose that the strategy of moni- 
toring RNA half-life will provide a powerful tool for investigating hitherto functionally uncharacterized regulatory RNAs. 

[Supplemental material is available for this article.] 



Whole transcriptome analyses using tiling microarrays (Bertone 
et al. 2004) and deep sequencing (Nagalakshmi et al. 2008) have 
revealed huge numbers of novel transcripts, including long and 
short noncoding RNAs (ncRNAs). The ratio of noncoding to pro- 
tein-coding genomic regions increases as a function of develop- 
mental complexity (Mattick 2004), suggesting that revealing the 
functions of ncRNAs transcribed from noncoding genomic regions 
is important for understanding genome function in higher or- 
ganisms. The ncRNAs can be roughly classified into two groups: 
small transcripts, such as microRNAs and piwi-interacting RNAs 
(piRNAs), and long transcripts (Prasanth and Spector 2007). Al- 
though the biological importance of small ncRNAs has been 
documented in recent years, the physiological functions of long 
ncRNAs (IncRNAs) are poorly understood. Recently, significant 
efforts have been applied to reveal the function of IncRNAs. Several 
approaches have succeeded in identifying dozens of functional 
IncRNAs (Guttman et al. 2009). However, the biological functions 
of the vast majority of IncRNAs remain unclear. Thus, novel prop- 
erties that can distinguish functional ncRNAs from transcriptional 
noise are required. 

Numerous studies of mRNAs have revealed that changing the 
abundance of transcripts by regulated RNA degradation is a critical 
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step in the control of various biological pathways (Keene 2010). 
It has been estimated that the mRNA abundance of 5%-10% 
of human genes is controlled through the regulation of RNA sta- 
bility (Bolognani and Perrone-Bizzozero 2008). It has been pro- 
posed that the specific half-life of each mRNA is closely related 
to its physiological function (Lam et al. 2001; Yang et al. 2003; 
Raghavan and Bohjanen 2004; Sharova et al. 2009; Rabani et al. 
2011; Schwanhausser et al. 2011). Although mRNAs of most house- 
keeping genes have long half -lives, mRNAs of many regulatory genes, 
which encode proteins that are required for only a limited time in the 
cell — such as cell cycle regulators, factors responsible for responses to 
external stimuli, and regulators of growth or differentiation — often 
have short half-lives. Moreover, most transcriptionally inducible 
genes are disproportionately classified into the group of genes with 
rapid mRNA turnover. It is possible, therefore, that the RNA sta- 
bility of noncoding transcripts also reflects their functions. 

Traditionally, RNA decay has been assessed by blocking global 
transcription with transcriptional inhibitors, e.g., actinomycin D 
(ActD), and subsequently monitoring ongoing RNA decay over 
time. However, inhibitor-mediated global transcriptional arrest 
has a profoundly disruptive impact on cellular physiology and 
interferes with the precise determination of the RNA degradation 
rate (Blattner et al. 2000; Friedel et al. 2009). Here, we present 
a novel inhibitor-free method (5'-bromo-uridine immunoprecipi- 
tation chase, BRIC) that enables measurement of RNA decay un- 
der nondisruptive conditions. Determination of the half- lives of 
whole transcripts by BRIC, combined with multifaceted deep se- 
quencing (BRIC-seq), suggest that there is a relationship between 



22:947-956 © 2012, Published by Cold Spring Harbor Laboratory Press; ISSN 1088-9051/12; www.genome.org 



Genome Research 947 

www.genome.org 



Tani et al. 



the stability of ncRNAs, as well as mRNAs, and their physiological 
functions. 

Results 

Genome-wide determination of RNA half-lives by BRIC-seq 

The BRIC method labels endogenous transcripts with 5'-bromo- 
uridine (BrU) by its addition to cell culture media (Fig. 1A). Total 
RNAs containing BrU-labeled RNAs (BrU-RNAs) are then isolated 
from cells at sequential time points after removal of BrU from the 
culture medium, and BrU-RNAs are recovered by immunopu- 
rification (Fig. 1A; Supplemental Fig. S1A,B). The recovery ratios of 
BrU-RNAs to nonlabeled RNAs were 10-100:1 when cells were la- 
beled with 150 fxM BrU, indicating that most RNAs were effectively 
labeled and were suitable for assaying decreasing levels (Supple- 
mental Fig. SI A). We also confirmed that BrU-RNAs were recovered 
within a linear range under the experimental conditions used 
(Supplemental Fig. SIB). The isolated BrU-RNAs can be used as 
templates for reverse transcriptase-mediated cDNA synthesis for 
subsequent quantification of RNA, because the BrU-RNA template 
does not cause misincorporation by reverse transcriptase (Core 
et al. 2008). Importantly, BrU did not cause harmful effects com- 
pared with other uridine analogs, such as 4-thiouridine (Friedel 
et al. 2009) and 5-ethynyluridine (Fig. IB; Supplemental Fig. SIC; 
Jao and Salic 2008). Moreover, BrU did not alter RNA localization, 
in contrast to ActD, which significantly affects the nuclear speckle 
localization of MALAT1, a nuclear IncRNA (Fig. 1C). Thus, BrU is 
a suitable agent to label RNAs for monitoring RNA degradation 



under physiological conditions. To test whether the BRIC method 
can determine RNA decay under nondisruptive conditions, we 
compared the RNA decay rates of MALAT1 (Tano et al. 2010) in 
HeLa Tet-off (TO) cells using the BRIC method and the Tet-off 
system, which does not use transcription inhibitors. The half-life 
of MALAT1 measured by the BRIC method and reverse transcrip- 
tion-quantitative real-time polymerase chain reaction (RT-qPCR) 
(ti/2 = ~7.6 h) was similar to that measured by the Tet-off system 
and RT-qPCR (t 1/2 = ~7.0 h) (Fig. ID). In contrast, ActD inhibited 
the degradation of MAL ATI (t 1/2 > 12 h) (Supplemental Fig. S2). 
Thus, the BRIC method allowed us to measure RNA decay under 
physiologically undisturbed conditions in vivo. 

We investigated the half-lives of mRNAs and ncRNAs to val- 
idate the assumption that RNA half-life is associated with the 
function of the RNA. We analyzed more than 5 million RNA-seq 
tags obtained from BRIC-deep sequencing analysis (BRIC-seq) at 
several time points after BrU pulse-labeling of HeLa TO cells. We 
calculated genie representations using "reads per kilobase of exon 
model per million mapped reads" (RPKM) as a means of normal- 
izing for gene length and depth of sequencing (Mortazavi et al. 
2008). We then mapped the sequencing tags to the human refer- 
ence genome sequence (hgl8) using Eland as a mapping software, 
allowing no mismatch settings. We used several public databases 
for annotation (see Supplemental Tables S1-S7; Methods). Com- 
parison of RNA-seq tags from two independent BRIC-seq experi- 
ments showed high reproducibility for transcripts with an RPKM 
>1 (Supplemental Fig. S3). Transcripts were selected for analysis 
that had an RPKM value >1 at time 0 h and determined the half- 




Figure 1. 5'-Bromo-uridine immunoprecipitation chase (BRIC) for measuring RNA decay in physiologically undisturbed conditions in vivo. (A) 
Overview of BRIC. (B) A549 cells were treated with various uridine analogs at the concentrations indicated for 48 h. The number of viable cells was counted 
using the Cell Counting Kit-8 (Wako). The relative abundance (1 00% in untreated cells) is shown in the graph. Values represent mean ± SD obtained from 
four replicate experiments. (C) A549 cells were treated with 5'-bromo-uridine (BrU) or Actinomycin D (ActD) at the concentrations indicated for 7 h. RNA- 
FISH was used to detect MALAT1 (green). Nuclear speckle proteins, SRSF2 (red) and SRRM2 (magenta), were assayed by immunostaining. Cells were 
counterstained with DAPI (scale bars, 1 0 ixm). (D) Decay rates of MALAT1 were determined by BRIC and RT-qPCR (gray line) or by the Tet-off system and 
RT-qPCR (black line) in HeLa TO cells. Relative quantitative values at time 0 h were arbitrarily adjusted to 1 00%. Values represent mean ± errors obtained 
from duplicate experiments. The half-lives of MALAT1 using BRIC or the Tet-off system were 7.6 h and 7.0 h, respectively. 
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lives for 11,052 mRNAs (average half-life: -6.9 h) and for 1418 
annotated ncRNAs (average half-life: —7.0 h). The actual number 
of ncRNAs in our analysis was less than the number of ncRNAs 
assumed previously, because we selected highly and moderately 
expressed transcripts with an RPKM > 1. Thus, the analysis is biased 
toward highly and moderately expressed transcripts. Typical map- 
ping data of RNA-seq tags are shown in Figure 2. BRIC-seq allows us 
to monitor the degradation of the RNAs lacking a poly(A) tail, such 
as histone mRNAs (Supplemental Fig. S4A,B). Moreover, we could 
determine the stability of pre-miRNAs (Supplemental Fig. S4C,D). 
The relative distribution of ncRNA half-lives was similar to that of 
mRNA half -lives (Fig. 2D), suggesting that the stabilities of ncRNAs 
are regulated, as is the case with mRNAs. 

We then performed an analysis of Gene Ontology (GO) terms, 
the association of functional categories with differences in mRNA 
half-life distribution (Fig. 3; Table 1). GO categories associated with 
regulatory functions were significantly represented in the short- 
lived mRNA group (t 1/2 < 4 h). In contrast, GO categories associated 
with housekeeping functions were disproportionately represented 
in the long-lived mRNA group (t 1/2 ^ 4 h). These results indicate 
that mRNA half-lives are correlated with the functions of genes, 



consistent with previous reports (Lam et al. 2001; Yang et al. 2003; 
Raghavan and Bohjanen 2004; Sharova et al. 2009; Rabani et al. 
2011; Schwanhausser et al. 2011). We further analyzed the corre- 
lation between GO terms and mRNA half-lives without splitting 
them into two groups in terms of RNA half -life. We confirmed that 
several GO terms, such as regulation of transcription, are associ- 
ated with short-lived mRNA. In contrast, several GO terms, such as 
translation, are associated with long-lived mRNAs (Supplemental 
Fig. S5). We then investigated the relationship between ncRNA 
half-lives and functional categories (Table 2). We found that several 
ncRNAs (tRNAs, small nucleolar RNAs [snoRNAs], and small Cajal 
body-specific RNAs [SCARNAs]) Qady and Kiss 2001) that are in- 
volved in housekeeping functions have long half-lives (ti/ 2 ^ 4 h). 
In contrast, well-known regulatory IncRNAs, such as CDKN2B-AS1 
(Yap et al. 2010; Kotake et al. 2011), HOTAIR (Gupta et al. 2010), 
TUG1 (Khalil et al. 2009; Yang et al. 201 1), and GAS5 (Kino et al. 2010) 
have short half -lives (t 1/2 < 4 h). We hypothesize that the half -lives of 
ncRNAs correlate with functional characteristics and that IncRNAs 
with regulatory functions may be enriched in the short-lived ncRNA 
group. We designated this novel class of ncRNAs as Short-Lived 
noncoding Transcripts (SLiTs), and we identified 795 SLiTs. 
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Figure 2. Typical BRIC-seq data. (A) Typical mapping data of RNA-seq tags obtained from BRIC-seq 
using a Genome Analyzer (lllumina). The chromosomal locus of AK091 718 is shown. (B) The relative 
level of AK091 71 8 RNA remaining was determined by BRIC through deep sequencing. The numbers of 
sequencing tags at the various time points are shown in A The relative quantitative values at time 0 h 
were arbitrarily adjusted to 100%. (C) The expression of AK091718 in HeLa TO cells assessed by 
Northern blot hybridization. (D) Distribution of RNA half-lives determined by BRIC-seq. BRIC-seq de- 
termined the half-lives of 11,052 mRNAs and 1418 ncRNAs expressed in HeLa TO cells. (Black circles) 
The fraction of mRNAs; (white circles) the fraction of ncRNAs. Transcripts with a half-life >24 h are 
included in the category ">24." 



Bioinformatic analysis of SLiTs 

To further characterize the SLiTs, we per- 
formed a series of bioinformatic analyses. 
Tissue profiling for SLiTs was performed, 
and the Z-score was calculated using the 
Adult Tissue TSSseq database (Supple- 
mental Fig. S6; Supplemental Table S8; 
Yamashita et al. 2010, 2011). The Z-score 
averages of SLiTs and long-lived ncRNAs 
were 0.75 and 0.74, respectively. The re- 
sults suggest that the tissue specificities of 
SLiTs are almost the same as those of long- 
lived ncRNAs. We also examined whether 
the half-life of an ncRNA correlates with 
other factors as follows. (1) Evolutionary 
sequence conservation was examined 
using the public database from http:// 
hgdownload.cse.ucsc.edu/goldenPath/ 
hgl8/phastCons44way /vertebrate/ (Sup- 
plemental Fig. S7; Supplemental Table S9). 
The average phastCons scores of SLiTs and 
long-lived ncRNAs were 0.08 and 0.29, 
respectively, and the P-value was <0.001. 
The results suggest that the evolutionary 
sequence conservations of SLiTs are lower 
than those of long-lived ncRNAs. The low 
conservation of SLiTs may be the result of 
recent and rapid adaptive selection. (2) 
The absolute expression levels between 
SLiTs and long-lived ncRNAs were exam- 
ined. The averages of the absolute ex- 
pression levels of SLiTs and long-lived 
ncRNAs were 42.8 and 544.7 RPKMs, re- 
spectively, and the P-value was below 
0.001. These results indicate that the ex- 
pression levels of SLiTs were lower than 
those of long-lived ncRNAs. This is rea- 
sonable, because the SLiTs are expected to 
show low abundance because of their fast 
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Figure 3. Representative data of mRNA stabilities. Relative levels of mRNA remaining were determined by BRIC-seq. The relative quantitative values at 
time 0 h were arbitrarily adjusted to 1 00%. (A) GO categories associated with regulatory functions. (Left) ATF1 is involved in transcription factor activity; 
(center) PLXNC1 is involved in multicellular organismal development; and (right) GPR157 is involved in G-protein-coupled receptor protein signaling 
pathway. (B) GO categories associated with housekeeping functions. (Left) SNRPB is involved with the spliceosome; (center) TIM Ml 7B is involved with the 
mitochondrion; and (right) EIF3F\s involved in translation. 



degradation. (3) The frequencies of the AU-rich element (ARE), 
which is the most well-known determinant of mRNA stability and 
whose consensus core sequence is AUUUA ; were examined in SLiTs 
and in long-lived ncRNAs (Supplemental Table S10). The fre- 
quencies of AREs in SLiTs and long-lived ncRNAs were 86% and 
77%, respectively, and the P-value was below 0.001. The frequency 
of AREs in SLiTs is potentially higher than that in long-lived 
ncRNAs. Thus, we found unique characteristics of SLiTs compared 
with long-lived ncRNAs. These unique characteristics may reflect 
the regulatory functions of SLiTs. 

Identification of novel IncRNAs involved in cell proliferation 
among SLiTs 

We investigated the function of a number of SLiTs in an attempt to 
identify important regulatory ncRNAs based on RNA half-life. We 



first selected 36 uncharacterized SLiTs (t 1/2 < 4 h) that did not 
overlap regions containing known protein-coding or microRNA- 
containing genes. Some of these SLiTs might have been nonsense- 
mediated decay (NMD) targets (Mendell et al. 2004), because NMD 
is a mechanism for eliminating aberrant mRNAs harboring pre- 
mature termination codons and other genomic noise, such as in- 
active transposons and pseudogenes (He et al. 2003; Baker and 
Parker 2004) that are likely to be nonfunctional. Put simply, 
transcripts that are degraded by NMD have a low potential to be 
bona fide functional ncRNAs. Based on this hypothesis, we ex- 
cluded NMD-targeted transcripts from later analysis. To identify 
transcripts that were degraded through NMD, we quantified SLiT 
expression in cells with an abrogated NMD pathway by in- 
troducing a UPF1 siRNA (Supplemental Fig. S8A; Supplemental 
Table SI 1) and judged that a given transcript was targeted by NMD 
when its level was increased by >200% in the UPF1 knockdown 



Table 1. GO categories significantly associated with mRNAs of different stabilities 



Gene with short-lived mRNA 



Gene with long-lived mRNA 



GO term 



Definition 



Number of GO term 



P-value 



Number of GO term 



P-value 



0003700 Transcription factor activity 

0007275 Multicellular organismal development 

0007186 G-protein-coupled receptor protein signaling pathway 

0005681 Spliceosome 

0005739 Mitochondrion 

0006412 Translation 



511 5.1 x 10~ 20 1 84 1 

407 2.9 x 10 6 209 1 

155 2.8 x 10 3 79 1 

36 1 74 1.7 x 10 7 
297 1 454 8.2 x 10 24 

37 1 91 6.8x10 11 



The degree of enrichment for a given cluster and category was analyzed by the hypergeometric distribution (P-value). Total numbers of genes with short- 
lived and long-lived mRNAs were 5557 and 5495, respectively. The degree of enrichment for a given cluster and category was analyzed by the hyper- 
geometric distribution (P-value). 
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Table 2. Relationship between functions and stabilities of ncRNAs 
Short-lived ncRNA Long-lived ncRNA 



Definition Number of RNAs P-value Number of RNAs P-value 



rRNA 


0 


1 


4 


3.3 x 10~ 2 


tRNA 


7 


1 


16 


8.1 x 10 3 


snoRNA 


70 


1 


95 


3.3 x 10" 5 


SCARNA 


4 


1 


16 


7.2 x 10~ 4 



The degree of enrichment for a given cluster and category was analyzed 
by the hypergeometric distribution (P- value). Total numbers of short-lived 
and long-lived ncRNAs are 795 and 594, respectively. 



cells (levels of the endogenous NMD-targeted ncRNA, SNHG1 
[Smith and Steitz 1998], were increased by 258%). Ten of the 36 
transcripts were judged to be NMD-targeted transcripts and were 
excluded from further analysis. We then determined the cellular 
localization of the remaining SLiTs. Among these 26 transcripts, 
15, two, and nine transcripts were localized to the nucleus, cyto- 



plasm, or both, respectively (Supplemental Fig. S8B; Supplemental 
Table SI 2). RNAs retained in the nucleus are the least likely to be 
translated; therefore, we judged that they were likely to be bona fide 
ncRNAs. Transcripts localized in the cytoplasm were enriched in the 
polysome fraction (data not shown) and were excluded. As a result, 
we identified 15 uncharacterized SLiTs that were likely to be bona 
fide ncRNAs (Supplemental Table SI 3). These 15 SLiTs are longer 
than 200 nt, fulfilling an established criterion for IncRNA classifi- 
cation. The degradation and abundance of three SLiTs — BX648321, 
NR_024586, and BX537481— were altered by retinoic acid, a typi- 
cal physiological stimulant involved in cellular proliferation and 
differentiation (Fig. 4; Ross et al. 2000). The abundance of all three 
SLiTs increased to >150% after retinoic acid treatment. Inter- 
estingly, their half-lives were very short (t 1/2 < 2 h). ActD, which 
influences the stability of mRNA, also affects 76% of SLiTs (Sup- 
plemental Figs. S9, S10; Supplemental Table S14). These results 
suggest that the stability of SLiTs is tightly regulated and supports 
the notion that short half-life indicates tight regulation and the 
possession of a particular biological function. 





1 -I 1 . 1 . — 1 -I 1 . 1 . — 1 -I 1 1 1 r- 

02468 02468 02468 
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Figure 4. Degradation of several SLiTs was altered by all-trans retinoic acid (ATRA) in HeLa TO cells. (A) Altered abundance of SLiTs by ATRA except 
BC01 8860, whose level was below the detection limit. HeLa TO cells were untreated (gray bar) or treated with 1 0 mM ATRA (black bar). Values represent 
mean ± errors obtained from duplicate experiments. (B) Decay rates of 5 SLiTs, whose levels were increased >1 35% in ATRA-treated cells compared with 
control cells in A and of GAPDHwere determined by BRIC and RT-qPCR in control cells (solid circle and black bar) and in ATRA-treated cells (open circle and 
gray bar). The relative quantitative values at time 0 h were arbitrarily adjusted to 1 00%. Values represent mean ± SD obtained from triplicate experiments 
[(**) P< 0.01; (*) P<0.05, Student's t-test]. 
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Functions of the 15 SLiTs were assessed by depleting each 
transcript using siRNAs and determining the number of viable cells 
(Fig. 5 A). We could only significantly reduce the levels of five 
transcripts (Supplemental Table SI 3), because siRNA tends to be 
less effective on nuclear RNAs compared with cytoplasmic RNAs 
(Ideue et al. 2009). Cell growth was significantly reduced by de- 
pletion of NR_015389, AK091718, and AK055657 with two dif- 
ferent siRNAs (Fig. 5 A). The knockdown experiment with two 
different siRNAs can minimize the possibility of the off-targeting 
on the phenotype. Consistent with the expectation from the 
knockdown experiment, cell growth was significantly increased by 
overexpression of NR_015389, AK091718, and AK055657. These 
results show that the three SLiTs are involved in cell growth. Im- 
portantly no floating cells (dead cells) were seen by depletion of 
NR_015389, AK091718, and AK055657 under microscopic obser- 
vation (data not shown), suggesting that these SLiTs are modifying 
factors but are not essential for cell proliferation. Biochemical 
fractionation indicated that NR_015389 ; AK091718, AK055657, 
AK123363, and NR_033895 were predominantly enriched in the 
heavy fractions of the sucrose density-gradient analysis, along with 
nuclear IncRNAs, such as NEAT1 (Clemson et al. 2009; Sasaki et al. 
2009; Sunwoo et al. 2009) and MALAT1, which are components of 
nuclear bodies (Fig. 5B; Supplemental Fig. Sll; Hutchinson et al. 
2007; Miyagawa et al. 2012). These five SLiTs may localize to spe- 
cific nuclear domains by forming a dense RNA-protein complex, 
similar to NEAT1 and MALAT1. 
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Figure 5. Functional analysis of SLiTs. (A) HeLa TO cells were treated with a control siRNA or with 
siRNAs targeting indicated SLiTs. The number of viable cells was counted using the Cell Counting Kit-8. 
Two siRNAs were used for NR_01 5389, AK091 71 8, and AK055657 to rule out the possibility of an off- 
target effect of siRNA (black or gray bars indicate the first or second siRNA, respectively). The relative 
abundance (1 00% in control cells) is shown in the graph. Values represent mean ± SD obtained in four 
replicate experiments [(**) P < 0.01, Student's t-test]. (B) HeLa TO cells were treated with plasmid 
vectors as indicated. The number of viable cells was counted using the Cell Counting Kit-8. The relative 
abundance (1 00% with mock vector) is shown in the graph. Values represent mean ± SD in six replicate 
experiments [(***) p < 0.001, Student's t-test]. (C) The relative amount of RNA was quantified by RT- 
qPCR from nuclear fractions obtained by sucrose gradient centrifugation. 



Discussion 

In this study, we report a novel method for determining the half- 
life of an entire transcriptome. Moreover, by analyzing the re- 
lationship between RNA half-lives and functional category, we 
found that the half-lives of ncRNAs correlate with functional 
characteristics. We propose a new category of ncRNAs that have 
short half-lives. We named this novel class of ncRNAs, Short-Lived 
noncoding Transcripts (SLiTs). To the best of our knowledge, this is 
the first report to propose a new category of ncRNAs according to 
their half-life. Bioinformatic analyses revealed the unique charac- 
teristics of SLiTs compared with long-lived ncRNAs. We suspect 
that the unique characteristics of SLiTs may reflect their regulatory 
functions. Finally, using SLiT criteria, we identified three novel 
nuclear ncRNAs that are involved in cell proliferation and three 
novel nuclear IncRNAs with stabilities regulated by retinoic acid. 

The nucleoside analog, 4-thiouridine (4sU), was previously 
used to label endogenous RNAs and mRNAs, and half-lives were 
indirectly estimated by comparing the preexisting mRNAs and 
purified 4sU-containing mRNAs (Friedel et al. 2009; Rabani et al. 
2011; Schwanhausser et al. 2011). According to these reports, the 
median half-life of mRNAs in mammalian cells varied (40 min to 
~9 h). This discrepancy may be due to inherent limitations in 
indirect estimation of RNA degradation rates by comparing puri- 
fied 4sU-containing mRNAs with preexisting mRNAs. Indeed, the 
yield of purification of 4sU-containing mRNAs would dramatically 
alter the estimation. In this regard, directly 
chasing decreasing levels of labeled RNA 
would be more suitable for determining 
RNA turnover. 

SLiTs are different from previously 
identified groups of ncRNAs, such as TSSs 
(transcription start site-associated RNAs) 
(Seila et al. 2008), PASR (promoter-asso- 
ciated short RNAs), PALR (promoter-as- 
sociated long RNAs), or TASR (termini- 
associated short RNAs) (Affymetrix/Cold 
Spring Harbor Laboratory ENCODE Tran- 
scriptome Project 2009), because SLiTs 
do not flank active promoters or termi- 
nators of genes. SLiTs are clearly detect- 
able without depleting the nuclear RNA 
degradation pathway; therefore, they 
are different from CUT (cryptic unstable 
transcripts) and PROMPT (promoter up- 
stream transcripts) (Wyers et al. 2005; 
Preker et al. 2008), which are only de- 
tectable in the absence of the nuclear 
RNA degradation pathway. The majority 
of SLiTs differ from lincRNAs defined in 
previous reports (Guttman et al. 2009; 
Khalil et al. 2009), because the expression 
levels of most lincRNAs are very low, 
whereas the expression levels of most 
SLiTs are significant, perhaps reflecting 
their importance in biological processes. 

In this study, several nuclear SLiTs 
were silenced by siRNAs. Indeed, previous 
reports stated that nuclear-localized 
ncRNA such as RN7SK (Robb et al. 2005; 
Berezhna et al. 2006), MALAT1 (Tano 
et al. 2010; Miyagawa et al. 2012), and 
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NEAT1 (Clemson et al. 2009) could be silenced by siRNAs. There 
are at least two possibilities for silencing nuclear-localized RNA by 
siRNA: (1) The component of siRNA-induced silencing localized in 
the nucleus can function to cleave endogenous nuclear target 
RNAs in the nucleus (Robb et al. 2005). (2) The cytoplasmic RISC 
cleaves endogenous nuclear target RNAs at the M phase of the cell 
cycle when the nuclear envelope disassembles and nuclear com- 
ponents, including nuclear ncRNAs, are contacted by cytoplasmic 
components such as the RISC complex. 

We propose that SLiTs are a new classification of regulatory 
IncRNAs, and we suggest that ncRNA half-life is a key parameter for 
regulating the function of IncRNAs and for identifying uncharac- 
terized functional noncoding transcripts. 

Methods 

Cell culture and drug treatments 

A549, HeLa TO (Clontech), and HEK293T cells were grown in 
Dulbecco's modified Eagle's medium (DMEM) supplemented with 
10% fetal bovine serum (FBS) and antibiotics at 37°C in a humidified 
incubator with 5% C0 2 . Cells were treated with uridine (Sigma- 
Aldrich), 5'-bromo-uridine (BrU; Wako), 4-thiouridine (4sU; Sigma- 
Aldrich), or 5-ethnyluridine (EU; Invitrogen) for 48 h, and the 
number of viable cells in a 96-well plate were counted using a Cell 
Counting Kit-8 (Dojindo), in accordance with the manufacturer's 
instructions. 

Plasmid constructs 

To construct pGEM-Luc2, the cDNA encoding firefly luciferase 2 
(Luc2) was amplified using the primer sets listed in Supplemental 
Table S15 and pGL 4.13 SV40 as a template. The cDNA was then 
cloned into pGEM-T (Promega), according to the manufacturer's 
instructions. To construct the tetracycline-dependent transcrip- 
tional repression system oiMALATl, the cDNA encoding human 
MALAT1 was amplified using a human cDNA library as a template. 
The cDNA was then cloned into the pGEM-T vector (Promega), as 
described by the manufacturer. The resulting cloned human 
MALAT1 cDNA spanned nucleotides 493-8443, corresponding to 
the reference sequence (RefSeq) accession number NR_002819. 
The tet response element (TRE) spliced to firefly luciferase 2 (Luc2) 
cDNA was amplified using the pTRE-Luc2 vector as a template. The 
pTRE-Luc2 vector was constructed using the pTRE promoter re- 
gion, which was amplified using a pTRE vector (Clontech) as the 
template. This fragment was inserted at blunted EcoRV and StuI 
sites, located upstream of the luciferase cDNA in pGL 4.13 SV40 
(Promega). The amplified TRE+Luc2 fragment was digested with 
Apal and ligated at the Apal site, upstream oiMALATl, in pGEM- 
T-MALAT1. The resultant plasmid (pTRE-Luc2-MALATl) expressed 
MALAT1 spliced to Luc2 cDNA at its 5' end (Luc2-MALAT-1) whose 
transcription was repressed by tetracycline or doxycycline. To 
construct the expression vectors for NR_0 15389 and AK091718, 
the cDNAs encoding human NR_015389 and AK091718 were 
amplified using a human cDNA library as a template. The cDNAs 
were cloned into pcDNA3.1/hygro(+) (Invitrogen), according to the 
manufacturer's instructions. The AK055057/FLJ31095/IMR321000193 
cDNA clone was purchased from NBRC, NITE. Additional infor- 
mation can be provided upon request. 

BRIC 

For pulse-labeling of newly transcribed RNA with BrU, typically, 
cells were incubated for 24 h at 37°C in the presence of 150 |jlM BrU 
in a humidified incubator with 5% C0 2 . At indicated time points 



after the removal of BrU-containing medium, cells were harvested 
for RNA preparation using RNAiso Plus (TaKaRa), and total RNA 
was isolated. 

To normalize the isolation of BrU-labeled transcripts between 
time points, BrU-labeled luciferase RNA was used as a spike-in 
control to serve as an internal standard. The BrU-labeled luciferase 
RNA was transcribed with 0.3 mM BrUTP, 7.5 mM rNTPs, and 1.6 
juLg of linearized p-GEM Luc2 using the T7 RiboMax Express Large 
Scale RNA Production System (Promega). The BrU-labeled luciferase 
RNA was purified using the QIAGEN RNeasy MinElute Cleanup Kit 
(QIAGEN), in accordance with the manufacturer's instructions. 

Anti-BrdU monoclonal antibody (mAb)-coated protein 
G-Sepharose beads were prepared as follows. Protein G-Sepharose 4 
Fast Flow (GE Healthcare) slurry (20 (jlL per sample) was washed 
three times with 1 mL of 0.1% bovine serum albumin (BSA) in 
phosphate-buffered saline (PBS). The final pellet was resuspended 
in 300 |xL of 0.1% BSA in PBS, and 1 fig of the anti-BrdU mAb 
(MBL), which also recognizes BrU, was added per sample. The 
mixture was incubated with rotation for 1 h at room temperature. 
The beads were washed three times with 1 mL of 0.1% BSA in PBS. 
The final pellet was resuspended in 200 jjlL of 0.1% BSA in PBS 
containing 100 U of RNasin plus RNase inhibitor (Promega). 

Two micrograms (for RT-qPCR) or 12 |xg (for deep sequencing) 
of BrU-labeled total RNA, or 0.8 ng of BrU-labeled luciferase RNA as 
a spike-in control, were prepared as described above, and dissolved 
in buffer A (25 mM Tris-HCl at pH 7.4, 6.25 mM EDTA). Each 
mixture was denatured by heating for 1 min at 80°C and then 
added to the anti-BrdU mAb-conjugated beads. The mixture was 
incubated for 1 h at room temperature with rotation. Beads were 
washed three times with 1 mL of 0.1% BSA in PBS. The final pellet 
was resuspended in 100 |xL of buffer B (10 mM Tris-HCl at pH 7.4, 
50 mM NaCl), and then 300 of ISOGEN LS (Nippon Gene) was 
added, followed by RNA isolation in accordance with the manu- 
facturer's instructions. The isolated RNA was used for RT-qPCR or 
deep sequencing. 

Deep sequencing and data analysis 

Using 1 juLg of RNA, an RNA-seq library was constructed using the 
mRNA Seq Sample Preparation Kit, according to the manufac- 
turer's instructions (Illumina). Thirty-six base-pair single-end-read 
RNA-seq tags were generated using an Illumina GA sequencer, 
according to the standard protocol. The fluorescent images were 
processed to sequences using the analysis Pipeline supplied by 
Illumina. RNA-seq tags were mapped to the human reference ge- 
nome sequence (hgl8) using Eland as a mapping software and 
allowing no mismatch settings. Splitting RNA-seq tags represent- 
ing exon-exon junctions were rescued according to Ref-seq ex- 
tracted from the UCSC Genome Browser. We calculated genie 
representations using "reads per kilobase of exon model per mil- 
lion mapped reads" (RPKM) as a means of normalizing for gene 
length and depth of sequencing (Mortazavi et al. 2008). RNA-seq 
tags were assigned to corresponding RefSeq transcripts or to 
IncRNAs when their genomic coordinates overlapped. We then 
analyzed two replicates of HeLa TO cells expressing BrU-incorporated 
RNAs using the Illumina high-throughput sequencing platform. 
The origin and the orientation of the RNAs were documented ge- 
nome-wide by mapping the reads to the reference human genome. 
In this study, we used criteria such that data >1 RPKM can be used 
to determine the half-lives of transcripts. We determined the half- 
lives of the RNAs by calculating the time when the RNA-seq value 
reached half of the initial RNA-seq value (0 h time point) (Sup- 
plemental Fig. SI 2). Data from time points after the time point at 
which the RNA-seq value reached half that of the initial value were 
omitted from the estimation. For instance, the 24-h time-point 
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value of the A RNA model and the 12- and 24-h time-point values 
of the B RNA model in Supplemental Figure S12 were omitted from 
our analysis procedure. By following this strategy, we could obtain 
the most accurate estimation, because the RNA-seq values of later 
time points with relatively inaccurate data were omitted from the 
calculation. To clarify the errors associated with RNA half -life mea- 
surements, we determined the error value, calculated (ave2 - avel)/ 
avel values for each transcript, and drew the distribution of these 
values in a graph of abundance frequency (Supplemental Fig. S13). 

Databases 

We used mRNA sequences available from public databases: 18,001 
RNAs from NM of RefSeq (Sayers et al. 2010). We selected used 
uncharacterized ncRNA candidates available from several public 
databases: 2829 RNAs from NR of RefSeq (Sayers et al. 2010), 1738 
RNAs from FLJ full-length human cDNAs (Ota et al. 2004; Wakamatsu 
et al. 2009), 36 RNAs from the Mammalian Gene Collection (MGC 
Project Team et al. 2009); Deutsches Krebsforschungszentrum 
(DKFZ); Kazusa DNA Research Institute (KIAA), and 261 RNAs from 
human large intergenic noncoding RNAs (Khalil et al. 2009). The 
sum of these RNAs was 5507, and 1418 out of 5507 RNAs passed 
the cutoff of RPKM > 1. We examined the long ncRNA databases, 
831 RNAs from human long noncoding RNAs Qm et al. 2010), and 
223 out of 83 1 RNAs passed the cutoff. We also examined the other 
IncRNA database, 5445 RNAs from the IncRNA database (http:// 
www.lncrnadb.org), and 1386 out of 5445 RNAs passed the cutoff. 
(See also Supplemental Tables S1-S7.) 

Reverse transcription-quantitative real-time polymerase chain 
reaction [RT-qPCR) 

The isolated RNA was reverse-transcribed into cDNA using the 
Quantitect Reverse Transcription Kit (QIAGEN) or the PrimeScript 
RT Master Mix (Perfect Real Time; TaKaRa). The cDNA was ampli- 
fied using the primer sets listed in Supplemental Table SI 7. Lucif- 
erase or GAPDH was used for normalization in BRIC or in other 
experiments, respectively. SYBR Premix Ex Taq II (Perfect Real 
Time; TaKaRa) was used in accordance with the manufacturer's 
instructions. Quantitative real-time reverse transcription PCR 
analysis was performed using a Thermal Cycler Dice Real Time 
System (TaKaRa). 

Functional analysis 

To identify functional groups significantly over-represented among 
short- or long-lived transcripts, we compared the overall distribu- 
tion of RNA half-lives against the distribution of RNA half-lives for 
specific functional categories. For this purpose, GO terms were 
associated to each gene using loc2GO at NCBI. Half-lives of the 
genes that were associated with a particular GO term were calcu- 
lated for each group. The statistical significance of the difference 
in the half-lives between different gene groups was evaluated by 
Wilcoxon's signed rank test. 

Tet-off system 

HeLa TO cells were transfected with plasmid using Lipofectamine 
2000 (Invitrogen), according to the manufacturer's instructions. 
Briefly, plasmid DNA was used at a final concentration of 1 |xg/mL. 
The cells were incubated for 24 h at 37°C in a humidified incubator 
with 5% C0 2 . After the addition of 1 |xg/mL doxycycline (tetra- 
cycline analog), cells were harvested at the indicated time points. 
Total RNA was isolated according to the method of Chomczynski 
and Sacchi (1987). To distinguish Tet-off promoter-driven exoge- 



nous MALAT1 from endogenous MALAT1, the amount of lucifer- 
ase cDNA (as a tag sequence) conjugated to the 5' end of MALAT1 
was measured by qPCR. 

RNA fluorescence in situ hybridization [FISH) 
and immunocytochemistry 

To prepare template cDNA for synthesis of the RNA probe, a 3.3-kb 
fragment of MALAT1 was amplified using the primer sets listed 
in Supplemental Table SI 7. The amplified MALAT1 fragment 
was cloned into pGEM-T (Promega). A digoxigenin (DIG)-labeled 
MALAT1 RNA probe was prepared using DIG RNA labeling mix 
(Roche) and T7 RNA polymerase (Promega), according to the 
manufacturers' instructions. The transcribed DIG-labeled RNA 
probes were partially digested in alkaline buffer (60 mM Na 2 C0 3 , 
40 mM NaHC0 3 at pH 10.2) (Cox et al. 1984). FISH was performed 
to visualize endogenous RNAs according to a previously described 
protocol (Sone et al. 2007). The hybridized probes were detected 
with standard immunocytochemical procedures, using the fol- 
lowing antibodies: anti-DIG sheep monoclonal antibody (Roche) 
and Alexa Fluor 488-conjugated anti-sheep IgG polyclonal goat 
antibody (Invitrogen). To visualize nuclear proteins, we used the 
following antibodies: anti-SRSF2 mouse mAb (Sigma- Aldrich) and 
anti-SRRM2 rabbit polyclonal antibody (kindly donated by Dr. 
Blencowe, University of Toronto). The secondary antibodies used 
were Alexa Fluor 594-conjugated anti-mouse IgG polyclonal goat 
antibody (Invitrogen) and Alexa Fluor 647-conjugated anti-rabbit 
IgG polyclonal goat antibody (Invitrogen). For the detection of 
DNA, samples were incubated with DAPI (Wako). Fluorescent im- 
ages were obtained using a confocal microscope TCS SP5 with 
a 63X/1.4 NA HC PL APO CS lens (Leica). Acquired confocal im- 
ages were pseudocolored and merged using the image analysis 
software LAS AF (Leica). All images shown are representative of 
several repeated experiments. 

siRNA treatments 

The sequences of the siRNAs used are listed in Supplemental Table 
SI 6. These siRNAs were transfected into cells using Lipofectamine 
RNAiMAX (Invitrogen), according to the manufacturer's instruc- 
tions. Briefly, siRNA duplexes were used at a final concentration of 
10 nM. When single transfection was performed, cells were har- 
vested 48 h or 72 h after transfection. When double transfections 
were performed, cells were retransfected 24 h after the first trans- 
fection, and harvested 72 h or 96 h after the first transfection. The 
qPCR analysis was used to determine whether RNA interference 
achieved significant depletion of each target sequence. 

Overexpression of SLiTs 

The expression vectors were transfected into cells using Lipofectamine 
2000 (Invitrogen), according to the manufacturer's instructions. 
Briefly, the expression vectors were used at a final concentration 
of 1 |xg/mL. The number of viable cells in a 96-well plate was 
counted 48 h after transfection using a Cell Counting Kit-8 
(Dojindo) in accordance with the manufacturer's instructions. 
Cells were also harvested 48 h after transfection. qPCR analysis 
was used to determine the overexpression levels of the SLiTs. 

Cell fractionation 

We followed the protocol of the Hirose Laboratory with slight 
modifications. HeLa TO cells (~1 X 10 7 cells) were collected using 
a cell scraper and centrifuged at SOOg for 5 min. The cell pellet was 
washed in ice-cold RSB150 buffer (10 mM Tris-HCl at pH 7.4, 150 mM 
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NaCl, 2.5 mM MgCl 2 ) and centrifuged. The cell pellet was resus- 
pended in 800 juuL of ice-cold RSB150 buffer. The cells were divided 
into two tubes for total or nuclear/ cytoplasmic fractions. For the 
total fraction, cells were resuspended in 0.5% Triton X-100 and 
placed on ice. For the nuclear/cytoplasmic fractions, 0.25 mg/mL 
digitonin was added to the cells and incubated for 5 min on ice. The 
cells were centrifuged at 3000^ for 1 min at 4°C, yielding the 
cytoplasmic (supernatant) and nuclear (pellet) fractions. The su- 
pernatant was reserved on ice as a source of the cytoplasmic frac- 
tion. The nuclear pellet was washed twice in ice-cold RSB150 buffer 
and centrifuged. The nuclear pellet was resuspended in 400 jjlL of 
ice-cold RSB150 buffer and 0.5% Triton X-100. The obtained 
fractions were subjected to RNA extraction using ISOGEN LS 
(Nippon Gene), in accordance with the manufacturer's instructions. 
Samples were eluted with 40 uX of RNase-free water for each frac- 
tion. The RNA solutions were treated with DNase I (TaKaRa); then 
a 10% RNA solution from each fraction was reverse-transcribed 
into cDNA and quantified by qPCR. 

Northern blot analysis 

Total RNA (10 u,g) was resolved using 1% formaldehyde agarose gel 
electrophoresis under denaturing conditions and then transferred 
to a nylon membrane. RNA blots were stained with methylene 
blue to check for equal loading and transfer. Hybridization was 
performed with ULTRAhyb (Ambion) using random primed 
[a- 32 P]dCTP-labeled cDNA probes. Hybridization signals were 
visualized and quantified with an FLA-9000 (FUJIFILM). 

Sucrose gradient centrifugation 

For cytosolic fractionation, HeLa TO cells (~5 X 10 6 per 10-cm 
dish) were washed twice with ice-cold PBS containing 100 |xg/mL 
cycloheximide and lysed directly on the plate by the addition 
of hypotonic lysis buffer (100 |xg/mL cycloheximide, 1 mM 
dithiothreitol, 200 U/mL RNase inhibitor [Promega], 1.5 mM KC1, 
2.5 mM MgCl 2 , 5 mM Tris-HCl at pH 7.4, 1% Triton X-100, and 1% 
deoxycholate). The lysate was collected and treated as described 
(Jefferies et al. 1994). Briefly, 500- uX aliquots were layered onto 
linear sucrose gradients (10%-50% in 80 mM NaCl, 5 mM MgCl 2 , 
20 mM Tris-HCl at pH 7.4, and 1 mM dithiothreitol) and centri- 
fuged at 50,000 rpm for 90 min at 4°C. The gradient was then 
fractionated using a Piston Gradient Fractionator (BioComp). RNA 
from each fraction was isolated using RNAiso Plus (TaKaRa) and 
analyzed by RT-qPCR, as described. For nuclear fractionation, iso- 
lated nuclei were disrupted by brief sonication, followed by cen- 
trifugation. The nuclear lysate was layered onto linear sucrose 
gradients (10%-50% in 80 mM NaCl, 5 mM MgCl 2 , 20 mM Tris- 
HCl at pH 7.4, 1 mM dithiothreitol) and centrifuged at 50,000 rpm 
for 90 min at 4°C. The fractionation of the sucrose gradient and 
isolation of RNA are described above. 

Data access 

The short-read sequence archive data that appear in this paper are 
registered in GenBank (http://www.ncbi.nlm.nih.gov/genbank) 
/DDBJ (http://ddbj.sakura.ne.jp). The data of two independent BRIC- 
seq experiments, with or without ActD, at three time points are 
registered as follows: the accession numbers DRA000357 (0 h, 
without ActD), DRA000358 (6 h, without ActD), DRA000360 (12 h, 
without ActD), DRA000361 (6 h, with ActD), and DRA000359 (12 h, 
with ActD). The data for determining the RNA half -lives at five time 
points are registered as follows: accession numbers DRA000345 (0 h), 
DRA000346 (4 h), DRA000347 (8 h), DRA000348 (12 h), and 
DRA000350 (24 h). 
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